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Abstract 

This paper describes a videogrammetric technique lor 
determining aerodynamic loads based on optical elastic 
deformation measurements. The data reduction methods 
are developed to extract the normal force and pitching 
moment from beam deformation data. The axial lorce is 
obtained by measuring the axial translational motion ol a 
movable shaft in a spring/bearing device. Proof-of-concept 
calibration experiments are conducted to assess the 
accuracy of this optical technique. 

1. Introduction 

Internal strain gauge balances have been used for years 
as a standard technique for measuring the integrated 
aerodynamic forces and moments in wind tunnels. A 
variety of internal strain gauge balances have been 
developed and the technical aspects of various balances 
have been studied in detail 1 1 1. Generally speaking, the 
structure of an internal strain gauge balance is complicated 
and the cost of fabrication is high. This paper presents an 
exploratory study for remotely measuring aerodynamic 
loads using a videogrammetic system. Unlike strain 
gauges, this method optically measures beam deformation 
to determine the normal force and pitching moment. The 
axial force is obtained by measuring the translational 
motion of a movable shaft in a spring/bearing device. 
Mathematical models for data reduction are developed to 
extract the aerodynamic forces and moments from 
deformation data. Uncertainty analysis is given to evaluate 
the contributions from the elemental error sources and 
correlation terms. At this stage, the normal force, pitching 
moment, and axial force are the primary quantities to be 
determined. In principle, the side force, rolling moment, 
and yawing moment can be determined in a similar manner. 
Proof-of-concept laboratory experiments have been 
conducted to validate the proposed methodology of 
measuring the aerodynamic loads. Potentially, the optical 
method can be used as an alternative to strain gauge 
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balances. In addition, the technique described in this paper 
can be easily integrated with the optical model attitude 
measurement techniques [2,3]. 


2. Cantilever Beam Deformation 

Deformation of a cantilever beam is utilized to 
calculate the normal force and pitching moment. Consider 
a cantilever beam bent by a force F and a moment M 
applied at the end, as shown in Fig. 1. According to the 
engineering beam theory in which the deformation due to 
shear strain is assumed to be negligible |4), the normal 
coordinate v of the beam is described by 

Eli x )v xx =-F( L-x )-M . (1) 

where v u . is the second derivative with respect to a, M is 
the local bending moment, E is the modulus of elasticity, 
and l(x) is the moment of inertia. The boundary conditions 
arc v v ( x -0)- i ' x0 and v( x = 0 ) - v„ . The solution for 

the displacement <S v - v - v n - v x() x is 


Sv< x) = -E~'( FL + M )Q t (x)+ E'FQJx). (2) 

where 

Q,(x)= f dx f / I(x" )I~ l dx" , 

Jo J o 

Q,(x)= f dx[ l I(x” )r i x"dx” . 

Jo Jo 

For a rectangular cross-section beam, the moment of inertia 
is l — wh * / 12 , where w is the width and h is the height of 
the cross-section. For a circular cross-section beam, the 
moment of inertia is / = nR 4 / 4 , where R is the radius of 
the circular cross-section. When a beam has a constant 
cross-section and the moment of inertia is constant, the 
expression (2) becomes 

6v(x) = -(Eir l (FL+M)x 2 /2+(EI)-' F x J /6 . (3) 

Therefore, the local displacement <5v and the slope change 


(4) 


Sv x =v ( -iv, are 

V / 6 El — x 2 / 2 El V F 
a - / 2EI -x/EI j{M t 
where M. =FL + M is the moment with respect to a 
moment center c. Eq. (4) gives a linear relation between 
the deformation (S\\8v x ) and the force and moment 
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(F,M C ). When the deformation due to shear strain is 

taken into account, a more complete analysis of a 
rectangular cantilever beam gives a non-linear relation [5] 

Sv = (x J /6EI )F-(x/2El )FL 2 +(I/3EI)FL < ‘ 

5v v =< x 2 / 2 El )F-( 1 / 2E1 )Fl} 

A generalized relation can be written as 

8v = ^C, j FL> and 8v x =^C 2j FL J . (6) 

j=o i=0 

where C # y are the coefficients to be determined in 

calibration tests using the least-squares method. These 
theoretical results provide basic models for data reduction 
to recover the force and moment. Generally, given the 
force and moment, the theoretical relations describe beam 
deformation profiles well. However, inversion to the force 
and moment from deformation measurements is a ‘stiff 
problem that is very sensitive to small errors. In practical 
data-reduction procedures, these theoretical results will be 
combined with empiricism to deal with imperfections in the 
real measurements. 

3. Data Reduction Methods 

3.1. Method based on local displacement and slope 

Eqs. (4) and (5) indicate that the force F and moment 
M ( depend on the local displacement Sv and the change 
of slope 8\\. Therefore, (F,M c ) can be determined 
from measurements of (8\\Sv x ). Here, the pitching 
moment M { is defined as M c - F( L~ L ( ), where L is the 
distance between the force (load) location and the beam 
support and L ( is the distance between the moment center c 
and the support. The local displacement and the change of 
slope are 

Sv = vi x ) - if x, vf ) - v x ( x rt , f )x rej 
8 i’ v = v x ( x)- v x ( x ref ) % (7) 

where x ref is the reference location where deformation is 
small such as the position near the support. Instead of 
directly using the definition (7), the practical procedures for 
calculating (5i\5v v ) are based on the affine 
transformation because of its robustness. 

The quantities Sv and 8v x are calculated from the 

optically measured coordinates of high-contrast targets on 
the beam. Figure 2 shows a typical layout of targets on a 
beam (sting) assembled to a wind tunnel model. A row of 
targets, denoted by B in Fig. 2, is placed in parallel to the 
beam centerline near the wind tunnel model. Another row 
of targets, reference targets denoted by A, is placed near the 
support of the beam. When the beam is deflected by 
aerodynamic loads, not only the targets at B move, but also 
the reference targets at A fnay move slightly since the 
support is not absolutely rigid. Thus, the total movement of 
target row B contains both the relative deformation at B 


with respect to A and the local movement at A. The local 
deformation at B, which is more sensitive to loads, can be 
obtained by removing the local movement at A from the 
total movement at B. 

In order to correct the movement of the reference target 
row A due to loads, the load-off (or wind-off) position of 
the target row A is used as a baseline position. Assuming 
that the local movement of the reference target row A can 
be approximated as rigid-body motion, one obtains the 
affine transformation between the load-off coordinate 


( X,Z ) and load-on coordinate ( X , Z ) for the reference 
target row A 
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where R A is the rotational matrix 


*A = 


cosO A sinO A 
—sinO A cosO A 


( 8 ) 


(9) 


0 A is a rotational angle at A, and T xA and T zA are the 
translations at A. The coordinate system ( X t Y ,Z ) is a 
standard wind-tunnel coordinate system in which X is in the 
free-stream flow direction, Y is in the spanwise direction of 
a wind-tunnel model at zero angle-of-attack, and Z is in the 
vertical direction. The beam deformation is in the ( X ,Z ) 


plane. The rotational angle 0 A and the translations (T u , 

T :A ) can be determined from the measured load -off and 
load-on coordinates of target row A by using the least- 
squares method. 

Applying the affine transformation (8) to the load-on 
coordinates of the target row B, we are able to eliminate the 
effect of the local movement at A. We denote the 
transformed load-on coordinates of target row B by 


( X Bon , Z Bon ) and call them the re-aligned load-on 

coordinates of target row B relative to the reference load- 
off target row A. The re-aligned load-on coordinates 

( X Bon , Z Bon ) of the target row B are related to the 


corresponding load-off coordinates ( X Bop ,Z Bog ) by the 
affine transformation 

( - \ ( ~ \ 


X Bon 
Z Bon 
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X Hoff 
Z Bop 
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where R r is the rotational matrix 


R b - 


cos6 b sinQ h 
-sinOa cosO 
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The change of the bending angle at B is 
Ar B = ( T xB , T :B ) is the displacement vector at B. 


GO) 


( 11 ) 

and 

The 
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values of 0 B and Ar B = (T xB *T :B ) can be determined by 
using the least-squares method and they are the average 
quantities of local deformation for target row B. Therefore, 
the change of the slope at B due to loads is simply 

8\\ - tan( 6 fi ). (12) 

The normal displacement Sv at B due to loads is 

Sv -n B • Ar /f , (13) 

where n B is the unit vector normal to the beam axis at B. 

In reality, the relationship between (F,M { ) and 
(8\\8\\ ) is more complicated than theoretical prediction 
by the engineering beam theory. The relations for 
calculating ( F,M ( ) are generally expressed as 

F = f f (8v,5v x ) and M t = f 2 (8\\Sv A ). (14) 

In practice, we do not tend to globally lit the whole set ol 
calibration data to obtain the complete functional relations. 
Instead, for a given data point ( 8v.8v x ) , we use a local 
2nd-order polynomial to interpolate a group of neighboring 
calibration data points, that is, 

F = B t 8 and M t = B 2 8 . (15) 

where S is the deformation vector defined as 

S = fSv.Sv ,( Sv ) 2 .SvSv .(Sv ) 2 f . 

A A A 

and B f and B 2 are the coefficients determined by 
calibration. 

A simple model approach can be also used to recover 
(F, M i ) from given (8v,8v x ). It is found that for a 
suitably chosen moment center L . the moment 
M - F(L-L C ) is simply a linear function of either <5v 
or 5v v . This is a reduced case in which Sv is 

proportional to <5^ . In this case, the moment M . is given 
by a simple proportional relation 

M c = F( L-L i )-a( 8v F (lb) 

where a is a proportional constant determined by 

calibration. The best moment center a* = L ( to achieve the 
linearity given by Eq. (16) can be determined by an 
optimization scheme. In addition, an empirical relation 
between F and Sv is 

Sv = f + p f ( L- L )+P 2 (L~L ) 2 IF. (17) 
where P 0 . fi n and fi 2 are constants determined in 

calibration. Eliminating (L-L c ) in Eqs. (16) and (17). 
one obtains a 2nd-order algebraic equation tor F 

P 0 F 2 +(P,M c - 8v)F + p 2 M ; = 0. (18) 

In principle. F and M t can be obtained from Eqs. (16) and 
(18) for a given 8 v . However, Eq. ( 18) has two real roots 

that are often close each other. It is not easy to choose the 
correct solution without a good initial guess of F. This is a 
shortcoming of the simple model approach. 


3.2. Method based on global beam deformation profile 

From Eq. (3), one knows that the displacement <Sv 

along the beam axis can be described by the theoretical 
relation 

<5 v( a ) = a a 2 + h a 1 , (19) 

where a is the coordinate along the beam axis Sv is 
obtained using Eq. (13), and the coefficients a and b are 
related to F and M { . However, the relation between {a. b) 
and ( F, M ) is not as simple as that given by Eq. (3). The 
empirical relations are symbolically expressed as 

F = fj( a , b ) and M ( = f 2 ( a . b ) . (20) 

For a given data point < a, b). a local 2nd-order polynomial 
fit to a group of neighboring calibration data points is used 
to recover ( F, M ( ) . 

As an alternative, we also use a simple model approach 
to calculate ( F , M ) . The empirical functional relations 


are 

a - l a 0 +a t ( L — L ( . ) + a 2 ( L- L F I F 
b = lb n +VL-L, ) + h 2 (L-L ) 2 IF. (21) 
where (a 0 ,ai.a 2 ) ^nd (h 0 ,b I ,b 2 ) arc determined in 
calibration. The solution to Eq. (2 1 ) is 

L-L =( -c, ± Fi -T77 )A 2f, ) 

F = a/[ a () + a j( L — L c ) + a 2 ( L - L ( F j . 
where c 0 - b 0 a/b — a 0 , c } -b t a/ b- a , 
c 2 - b 2 a / b - a 2 . 

3.3. Uncertainty 

For uncertainty analysis, consider the 
functional relations between ( F, M ( ) and ( <5v. Sv x ) 

F = f,( Sv,Sv x .p ) and M t = f 2 ( 8v,8v x ,q ) , (23) 

where p = / p } , ' Ps I an ^ ^ = / ci i ,q 2 * tm9£ lN I are 

the parameters characterizing the functional relations. 
Standard uncertainty analysis [6] gives the error 
propagation equations for the relative variances of F and 
M c . 


( 22 ) 

and 


general 


var(F) 
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X s - 

i = l 
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var(8\\ x ) 
i Sv x ) 2 
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+ S 7 A 


var(M ( ) _ r: var (^i I 
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var{ Sv } ( c 2 
f Sr F 


(24) 
varf Sv x ) 


{ S\ y F 
lvar( Sv )var(Sr x )j !/ ~ 


8\8v , 


(25) 


where the sensitivity coefficients are 
S, r =( Pj/F )( dF/dp, ), 
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S 2Vl =( q ,/M t . )(dM,/dq, ), 

S ISv =( Sv/F )( r)F/r)Sv ) , 

S lSv =< 5\' X /F MdF/ddv , ), 

S 2 sv =( Sv/M c )( dM,/dSv ).and 
S 2 *, = (Sv x ./M l )( dM, /dSv x ). 

The correlation coefficient R( C, • f , ) between the 
variables C , and C, is defined as 

R( C, .Qj ) = cov( CiCj )/l var< f, )var( Q } )J M , and the 

variance and covariance are var(£ i ) = < A£ t 2 > and 
covfCjC j ) ~ < d£id£ j > , where the notation < > denotes 

the statistical assemble average and A£ f = £,-<£,- > is 
the variation. There are two types of contributions to the 
total uncertainties in F and M c . The summation terms in the 
right-hand side of Eqs. (24) and (25) are the uncertainties in 
fitting the coefficients p and q in the data-reduction 
models. The other three terms in the right-hand side are the 
uncertainties in measurements of (8v,8v x ). The 
uncertainties in measurements of (8\\8v x ) are 
determined by the accuracy of the camera calibration (about 
1:60000 for the videogrammetric system used). In addition, 
an important error source that is not included in Eqs. (24) 
and (25) is a bias error in the data-reduction mathematical 
model itself since the model may fail to describe correctly 
actual measurement processes. 

The last terms in the right-hand side of Eqs. (24) and 
(25) are the correlation terms that are more intriguing since 
they could be positive or negative. The total uncertainties 
decrease when the correlation terms are negative and 
otherwise the total uncertainties increase. Based on the 
linear theoretical relation (4), the correlation between Sv 

and 8v x can be calculated, that is, 

< A( 8 v )A( 8v x )> = 

x* < ( af y > , „ r r , , . . , * 

— [x-2( L-L ( )][ x-3( L-L c )] 

12( El ) 2 

The correlation terms in Eqs. (23) and (24) are 
r 0 < A(8v)A( 8v x )> _ 

' ovov K 

^AEL^(x-2(L-L c )Jl x-3( L-L ( )l 
x"F 

c c <A(Sv)A(Sv l )>_ 

t** 2 *' Sv5v x 

— — — - J 7 ~ l x — 2( L — L )][x-3(L-L c )] 

6M‘ 

The correlation terms are negative 
2( L- L c ) < x < 3( L - L c ) and 0 < L c < L . Furthermore, 


(27) 

for 


at a fixed target location x = rL ( r < / ), the correlation 
terms are negative when the moment center jc - L c is in 
( 2 - r )L/ 2 < L c <( 3 - r )L/ 3 . This analysis indicates 
that the moment center can be suitably selected to reduce 
the total uncertainties in measurements. 

4. Videogrammetric System 

In this study, elastic deformation of a beam is 
measured using a videogrammetric system. Based on the 
principles of close-range photogrammetry, the 
videogrammetric system measures the coordinates of 
targets distributed along the beam from target centroids in 
digital images. Deformation of the beam is then calculated 
from the measured target coordinates. Figure 3 shows a 
schematic of a two-camera videogrammetric system for 
deformation measurements. Basic hardware of the system 
consists of two Hitachi KP-F1 CCD cameras with 15-mm 
lenses, a Dell PC with a Matrox Pulsar frame grabber 
board, light sources and retro-reflecting targets on the 
beam. Software includes programs for image acquisition, 
target-tracking/centroid calculation, and camera calibration. 
This system is able to provide three-dimensional 
coordinates ( X , Y , Z) in almost real time. The accuracy of 
the videogrammetric system used for this work is typically 
1:60000. After the coordinates of the targets on the beam 
are obtained, the deformation ( 8v , 5\\) can be calculated. 
The data-reduction programs for calculating the 
deformation, force and moment are written in Matlab. The 
technical aspects of the videogrammetric system have been 
described in detail by Liu et al. [7,8]. Comprehensive 
reviews of application of videogrammetry to wind tunnel 
testing have been given by Burner & Liu [2], and Liu el al. 
[ 8 ]. 

5, Experimental Results 
5.1. Simple beams 

Figure 4 shows a standard balance calibration 
apparatus and a two-camera videogrammetric system. 
Three simple beams have been calibrated to examine the 
proposed methodology of recovering the normal force and 
pitching moment (see Fig. 5). One is a brass beam with a 
0.75 inx0.75 in square cross-section and the rigidity of 

El - 350x10 * lb-in 2 . A more flexible stainless steel beam 
is actually a flat plate that has a rectangular cross-section of 

J.5inx0.25in and the rigidity of El =47x10* lb-in 2 . A 
more rigid stainless steel beam was also used, which has a 
circular cross-section with a diameter 0.75in and the 

rigidity of El =442x10 3 lb-in 2 . A number of retro- 
reflecting targets used for deformation measurements are 
placed along the centerline of a beam. 

Figures 6 and 7 show the local displacement and slope 
change ( 8v , 8v x ) as a function of the normal force for the 
brass beam at four different loading positions. It is found 
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that the dependence of both the local displacement and 
slope change on the pitching moment is linear when the 
moment center is suitably chosen. The linear relations are 
clearly shown in Fig. 8 for the brass beam, where the 
moment center determined by an optimization scheme is 
located at L c = 9.03 in. This linearity is utilized in the 
simple model method for data reduction. Figure 9 shows 
the errors in measurements of the normal force and pitching 
moment by using the method based on local displacement 
and slope change for the brass beam. The local 2nd-order 
polynomial and the simple model methods are used lor data 
interpolation. For the soft and stiff stainless steel beams. 
Figure 10 shows the measurement errors in the normal 
force and pitching moment. Typically, the relative errors in 
the normal force and pitching moment are within ±10% and 
±5%, respectively. 

The method based on global deformation profile is 
used to recover the normal force and pitching moment for 
the brass beam. Figure 1 1 shows measured deformation 
profiles of the brass beam for different loads at the loading 
position L — 28.069 in. These measured data can be well fit 

by the theoretical solution 8 v( x ) = a x 2 + b x ' . The errors 
in measurements of the normal force and pitching moment 
are shown in Fig. 12, where the local 2nd-order polynomial 
and simple model methods are used for data interpolation. 
Compared to the method based on the local displacement 
and slope change, the method based on the global 
deformation profile gives smaller relative errors in the 
normal force (within ±5%) and pitching moment (within 
±3%) for the brass beam. 

5.2. Sting-model combination 

A steel sting-model combination used in the Unitary 
Tunnel at NASA Langley was calibrated, as shown in Fig. 
13. The sling has a tapering circular cross-section with a 
linearly decreasing radius from 0.5975in. to 0. 197in. 
Accordingly, the rigidity decreases from 2800'X 10 to 
38x10' lb-in 2 and the mean rigidity of the sting is 

930x10 * lb-in 2 . In calibration tests, the maximum 
displacement of the sting is about 0.03in and the maximum 
change of the local bending angle is about 0.23 degrees. 
Two strain gauge bridges were also installed on the sting 
for measurement of the normal force and pitching moment, 
allowing a direct comparison between the strain gauge 
method and the optical method. 

During calibration tests, a number of dead weights 
from 2 lb to 10 lb were loaded at three different locations 
on the model to generate the required forces and moments. 
Three retro-reflecting targets were placed near the model 
for measuring local deformation and four other targets were 
placed near the support as a reference. The local 
deformation quantities (<5v, 8\\ ) were measured using the 
videogrammetric system. The method based on the local 
displacement and slope change was used to determine the 
normal force and pitching moment. Figure 14 shows the 


measurement errors in the normal force and pitching 
moment for this steel sting-model combination. The results 
indicate that the optical method based on remote 
deformation measurement is less accurate than the more 
sensitive strain gauges. The relative errors in the normal 
force and pitching moment obtained by the optical method 
are about ±5% in comparison with ±2% given by the strain 
gauges. 

5.3. Three-force beam 

In order to measure the axial force along with the 
normal force and pitching moment, a three-force beam has 
been designed and fabricated. As shown in Fig. 15. the 
three-force beam consists of a simple beam and a 
spring/bearing device that only allows translational motion 
along the beam axis. The structure of the spring/bearing 
device is shown in Fig. 16. A steel axial load shaft (made 
of medium-alloy A-2 steel) is a moving rod that is 
constrained by a linear bearing mounted inside the 
hardware assembly. The rod is hardened and has a ground 
surface finish of 0.625 (+().()()()2/-0. ()()()()) inches on the 
diameter. This tight tolerance allows a very close interface 
fit with the bearing so all lateral movement is minimized. 
Another main component of the device is the spring that 
balances the applied axial force. Ranges of the springs with 
different spring constants are available and can be selected 
to meet the requirements of the axial force calibration. The 
three-force beam attaches to the extended end ol the 0.625 
diameter stainless steel beam. The device housing has an 
outside diameter of 2.38 inches and is 6.0 inches long with 
an extended length of the entire beam assembly of 28.25 
inches in the balance calibration apparatus. As shown in 
Fig. 15, a number of retro-reflecting targets are placed on 
the simple beam for measuring beam deformation. Four 
targets are placed on the spring/bearing device as a 
reference and the other four targets are placed on the 
movable shaft for measuring the axial translational motion 
relative to the reference targets. The average spacing 
between retro-reflective targets is approximately 0.65 
inches. 

The normal force and pitching moment are obtained by 
using the same methods as previously described. The 
measurement errors in the normal force and pitching 
moment are shown in Fig. 17, indicating ±4% errors for the 
normal force and ±3% errors for the pitching moment. The 
axial force F (IX . is related to the axial displacement <5.v by a 
linear relation 

^ -F ilx n=kSx, (28) 

where k is an effective spring coefficient and F ax0 is the 
force at 8x = 0 . When the normal force F N and pitching 
moment M exist, k and F ax0 are not constants and they 
are weakly dependent on F N and M t . due to the presence 
of friction in the bearing. Experiments show that k and 
F ax0 are mainly related to the loading position M ( / F N 
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over a certain range of the normal force. Thus, an empirical 
calibration relation for the axial force is 

F IU =k( M, / F N )8x + F ilx0 (M c / F N ) , (29) 

where k( M ( / F N ) and M c / F n ) are empirically 

expressed by polynomials. Figure 18 shows the linear 
relation between the axial force and the measured axial 
displacement. As shown in Fig. 19, the measurement errors 
in the axial force obtained by using Eq. (29) are within 
± 8 %. 

6. Conclusions 

The methodology of optically measuring aerodynamic 
loads is developed based on beam deformation theory. A 
two-camera videogrammetric system is used for optical 
deformation measurements. The data reduction models for 
extracting the normal force and pitching moment utilize 
either the local displacement and slope change or the global 
beam deformation profile. The interpolation methods for 
calibration data include the local 2nd-order polynomial fit 
and the simple model approach. On the other hand, the 
axial force is obtained by measuring the translational 
motion of a movable shaft in a spring/bearing device. In 
order to validate the proposed technique, three simple 
metallic beams with different rigidities have been tested to 
recover the normal force and pitching moment. Typically, 
the relative errors in the normal force and pitching moment 
are, respectively, from ±5% to ±10% and ±2% to ±5% for 
the three beams. A steel sting-model combination was also 
calibrated for a direct comparison between the strain gauge 
method and the optical method. The loads obtained by the 
optical method have larger errors than the strain gauges. To 
determine the axial force, a three-force beam that allows the 
axial translational motion has been designed and fabricated. 
The measurement errors in the axial force are within ±8%. 
At this stage, the optical method is less accurate than the 
we 11 -developed strain gauge method. Nevertheless, the 
optical method has the capability of remote non-contact 
measurement and it can be readily integrated with the 
optical methods of model attitude measurement. The 
optical method described in this paper can be considered as 
an alternative in certain cases where the strain gauge 
method cannot be easily applied. Further research effort 
will be made to improve the accuracy of the optical method 
in measuring aerodynamic loads and apply the technique to 
wind tunnel testing. 
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Fig. 7. Local slope change as a function of the normal force 
at four loading positions for a brass beam. 
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Fig. 8. Local displacement and slope change as a function 
of the pitching moment for a brass beam. 
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Fig. 9. Relative errors in (a) the normal force and (b) 
pitching moment for a brass beam when the data reduction 
method is based on local displacement and slope change. 

0.3 


CD *=• 


0 

CL 


0.2 

0.1 

0.0 

- 0.1 

- 0.2 

-0.3 

-0.4 


0.15 

0.10 

0.05 

0.00 

-0.05 

- 0.10 

-0.15 

- 0.20 


• Local 2nd-order polynomial, soft steel beam 

* Simple model method, soft steel beam 

o Local 2nd-order polynomial, stiff steel beam 
a Simple model method, stiff steel beam 


10 


Force (lb) 

(a) (see caption in Fig. 10) 






.go 8 

** 8*0 

A A S A 


0 0 


O 0 0 0 

a • 

▲ o 


• Local 2nd-order polynomial, soft steel beam 

* Simple model method, soft steel beam 

o Local 2nd-order polynomial, stiff steel beam 
a Simple model method, stiff steel beam 


20 


40 


60 


80 


100 120 


160 


Moment (lb-in) 

(b) 

Fig. 10. Relative errors in (a) the normal force and (b> 
pitching moment for a soft steel beam and a stiff steel beam 
when the data-reduction method is based on the local 
displacement and slope change. 
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Fig. 1 1 . The measured deformation profiles of a brass beam 


and the fit by S if x ) = a x 2 + b x J at the loading position L 
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Fig. 1 2. Relative errors in (a) the normal force and (b) 
pitching moment for a brass beam when the data-reduction 
method is based on the global deformation profile. 



Fig. 13. The sting-model configuration with targets. 
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Fig. 14. Comparison between the optical method and strain 
gauges in measurements of (a) the normal force and (b) 
pitching moment for the steel sting-model configuration. 
The optical method is based on the local displacement and 
slope change. 
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Fig. 15. Three-force beam with targets. 



Fig. 16. Structure of the spring/bearing device. 
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Fig. 1 8. The axial force as a function of the axial 
displacement. 
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Fig. 19. Measurement errors in the axial force for the three- 
force beam at different normal forces and loading positions. 
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Fig. 17. Relative errors in (a) the normal force and (b) 
pitching moment for the three-force beam at different axial 
loads. The data-reduction method is based on the local 
displacement and slope change. 
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